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Abstract. We present a quantitative criterion for the thermal stability of thermonuclear shell sources. We find the thermal 
stability of shell sources to depend on exactly three factors: they are more stable when they are geometrically thicker, less 
degenerate and hotter. This confirms and unifies previously obtained results in terms of the geometry, temperature and density 
of the shell source, by a simplified but quantitative approach to the physics of shell nuclear burning. We present instability 
diagrams in the temperature-density plane for hydrogen and helium shell burning, which allow a simple evaluation of the 
stability conditions of such shell sources in stellar models. The performance of our stability criterion is demonstrated in various 
numerical models: in a 3 Mq AGB star, in helium accreting CO white dwarfs, in a helium white dwarf which is covered by a 
thin hydrogen envelope, and in a 1.0 Mq giant. 
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1. Introduction 

Thermonuclear fusion is the main source of energy in stars. 
To use fusion as main energy supply over long time scales re- 
quires extreme stability, since the thermonuclear reaction rates 
are sensitive functions of the temperature. I.e., a slight tempera- 
ture increase may enhance the nuclear burning rate drastically; 
if this results in an even higher temperature, a thermal runaway 
occurs. In most stars, such a runaway does not occur, since an 
increased energy output due to an enhancement of the nuclear 
burning rate leads to an overpressure and thus to an expansion 
of the burning region (Kippenhahn & Weigert 1990). However, 
two different circumstances have been shown to be able to pre- 
vent such a pressure increase. 

The first one is electron degeneracy. If nuclear burning 
starts in a degenerate gas, the increased energy production will 
lead to a higher temperature. But as — for complete degeneracy 
— the pressure does not depend on the temperature, it remains 
constant. This situation is realized during the core helium flash 
in low mass stars, at carbon ignition in Chandrasekhar-mass 
CO white dwarf, i.e. in a Type la supernova explosion, and in 
strong nova shell flashes. 

The second circumstance is a strong geometrical confine- 
ment of the burning region. As discovered by Schwarzshild & 
Harm (1965 1 and Weigert (1966 1, geometrically thin helium 
shell sources in red giants can become unstable. Since the lo- 
cal gravity in a thin shell source remains constant even when 
it expands in response to an increased energy production, hy- 
drostatic equilibrium will enforce a constant pressure in the 
shell. The thin shell instability is prone to occur in nuclear 



shell sources around compact cores, as the small scale heights 
induced by the compactness of the core ensures the geometri- 
cal confinement of the shells. Examples are thermal pulses in 
asymptotic giant branch stars (e.g. Iben & Renzini 1983 , Busso 
et al. 2001 1 Lugaro et al. 2003 1, hydrogen or helium shells on 



white dwarf s (e.g. Fujimoto fT982al[T 982b Cassisi et al. fl998l 
Langer et al. 2002 , Yoon et al. 2004 1, and X-ray bursts on neu- 
tron stars (e.g. Fujimoto et al. 1981 , Taam & Woosley 1996). 

In principle, both destabilizing circumstances may occur 
together. In central burning regions, however, the degeneracy 
aspect dominates, as gravity and thus pressure will not be con- 
stant in an expanding core. In nuclear shell sources, degen- 
eracy and geometrical confinement may simultaneously con- 
tribute to the instability of shell sources. Among various au- 
thors investigating the stability of nuclear shell sources (e.g. 
Giannone & WeigertfT967l Unno lT970l Den nis lT97ll Harm & 
Schwarzshild [19721 Stot hers & Chin 119721 Sackmann[l977j, 
Kippenhahn & Weigert J1990I ) developed a semi-quantitative 
understanding of the relevance of both destabilizing factors. 

The purpose of this paper is to develop the insight into the 
shell stability mechanism further such that a quantitative tool is 
obtained, which allows to assess the thermal stability of shell 
sources in numerical stellar models. This will allow for a phys- 
ical interpretation of fluctuating nuclear shell sources found 
in non-linear time-dependent stellar evolution calculations. As 
described in Section 2, we essentially follow the method by 
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Giannone & Weigert (1967), but assume homology in the re- 
gion of the shell source as in Kippenhahn & Weigert ( 1990). 
In Section 3, we derive thermodynamic parameters for partial 
degenerate conditions, which we use in Section 4 to formulate 
and investigate the stability criterion in idealized situations. In 
Section 5, we apply the criterion to various numerical stellar 
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models, and give a summary and discussion of our results in 
Section 6. 

2. Secular behavior of shell burning 

2. 1. Gravothermal specific heat 

Let us consider a shell source with a geometrical thickness D. 
The mass of the shell source is given as AM S = J B 4irr 2 pdr, 
where r is the radius of the bottom of the shell source and 
r s = r o + D, its upper boundary. Assuming r is constant, the 
relation between density and radius perturbation is given from 
dAAL = as 



S_p 

P 



3D/r s - 3(£>/r B ) 2 + {D/r s f r s 



(1) 



(Huang & Yu 1998). If we assume that the matter in the shell 
source expands or contracts homologously, we have SP/P = 
—4Sr/r from the hydrostatic equation and we obtain: 



SP 
~P 



Sp 



(2) 



where 



n s 



- (3D/r s - 3(D/r s ) 2 + {D/r s f) . (3) 

Kippenhahn & Weigert's value a s =4D/r s is restored when 
D/r s ^ 1, while it becomes 4/3 with D/r s =l, which is the 
case when the whole star contracts or expands homologously. 
Using the thermodynamic relation Sp/p = apSP/P — 
a T ST/T, it follows from Eq. (0 that 

SP a s ar ST Sp ax ST 

__ __ 



1 T 



(4) 



a s ap — 1 1 p a s ap 

where ctp = (d\np/d\nP)T and olt = — (d In p/d In T)p 
(Kippenhahn & Weigert 1990l. As we will see in the next sec- 
tion, ax and ap are functions of the degree of degeneracy and 
are not constant over a shell source. However, since we are in- 
terested in the mean properties of a shell source in this study, 
we will assume that like the density and the pressure perturba- 
tion, the temperature change is also homologous in the shell, 
and mean values over the shell source will be used for ap and 
ax- 

The heat perturbation in the shell source can be ob- 
tained from the first law of the thermodynamics as shown by 
Kippenhahn & Weigert | f!990t : 



Sq = c*ST 



with 



cp 1 - V 



a s a>T 



ad " 



(5) 



(6) 



a s ap — 1 

where cp denotes the specific heat at constant pressure and 
V ad = (dlnT/dlnP) ad . 

The quantity c* is called the gravothermal specific heat 
(Kippenhahn & Weigert 1990 1. When the heat Sq is added to 
the shell source, more energy than Sq is consumed by the ex- 
pansion work if c* < 0, thereby decreasing the internal energy. 
However, if c* > 0, the additional energy heats up the matter 
and thus favors instability. 



2.2. The stability criterion 

The gravothermal specific heat alone does not describe the ther- 
mal reaction of a shell source, as it is also influenced by its ther- 
mal interaction with its surroundings. The energy conservation 
in a star is expressed by the following equation: 



~dM r 



e N TT- 



dq 
dt' 



(7) 



For a shell source with the thickness D = r s — ro and the mass 
AM g , integration over the burning shell gives 



1j , 



L rQ — £n — -kg 



(8) 



L Tb and L ro are the luminosities at r — r s and r = ro, respec- 
tively. Ln is the luminosity due to the nuclear burning in the 
shell: Ln = J^m e ^>dM r where €n is the nuclear energy gen- 
eration rate in the shell source. L g is defined as J AM -^dM r . 
In AGB stars or accreting white dwarfs the main contribution 
to L g is due to the gravitational energy release by contraction. 
This term is not significant compared to Ln in general, and we 
will assume that the shell source is initially in a stationary state 
(i.e., dq/dt = 0). 

Perturbing Eq. (fJJ yields 



SL r = AALScn - AM, 



dSq 
~dt 



(9) 



where L ro ^C L Ts has been assumed. With the use of equa- 
tions (0 — 1|6) and with assuming radiative heat transfer, i.e., 



dT 



3 nL r 



dM r 64:n 2 ac r 4 T 3 ' 

we obtain the following equation for the temperature perturba- 
tion 9 = ST jT: 

T th 6 = ad, (10) 



with 



and 



Tth 



AM s Tcp 



v — 4 + kt 



OCT 



a s ap — 1 



(A + a s + k p ) 



c*/c p 



(11) 



(12) 



Here, kt = (dlnK,/dlnT) p , Kp = (d\nK/d\np)T, v = 
[d In e/d In T) p , and A = (d In e/d In p)r- The quantity Tth 
defined in Eq. (lilt corresponds to the thermal time scale of the 
shell source. Note that Eq. dlOt is essentially the same as the 



equation (16) in Giannone & Weigert ( 19671, for the case that 
the assumption of homology in the shell source is adopted in 
their analysis. 

When a > the shell source becomes thermally unstable 
and the perturbation growth time scale becomes 



Tgrowth = T"th/c. 



(13) 



The larger the value of a, the more rapidly the thermal insta- 
bility grows. If a < 0, nuclear burning in a shell source is ther- 
mally stable, and, upon a temperature increase, the temperature 
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drops faster for larger a. Therefore, a can serve as a measure 
of the susceptibility of a shell source to thermal instability. 

For the following discussions, let us define the numerator 
in Eq. d!2i as cte such that a — (Te/(c*/cp) and 



cte = v — 4 + kt 



OCT 



a s ap — 1 



(A + a s + k p ). 



(14) 



The first 3 terms in the above equation are related to the tem- 
perature perturbation and the remaining terms to the density 
perturbation. 

For a given c*, i.e. a given equation of state, the stability of 
a shell source is determined by the sign of cte. Physically, cte > 
means that the additional energy production in response to a 
positive temperature perturbation exceeds the additional energy 
loss by radiation. In other words, for ST > 0, we have Sq > if 
<te > 0, in which case a shell source is stable if c* < 0, since, 
the additional heat is consumed mostly for the expansion work 
of the shell source as discussed in section I27TI If cte < 0, the 
additional energy loss by radiation is larger than the additional 
energy production (i.e., Sq < 0) in response to ST > 0. If c* < 
and cte < 0, this results in further increase in the temperature 
perturbation since the net energy loss is compensated by the 
contraction of the system, increasing the internal energy of the 
shell source. However, this instability mode has little relevance 
as in stars cte hardly becomes negative if c* < 0, as we will 
see in Sect. 13.11 If c* > and cte < 0, the net energy loss 
results in the decrease in the internal energy and stable burning 
is established. 

In the following sections, we will investigate the possible 
modes of stability/instability of a shell source in a more quan- 
titative way. 

3. Physical conditions for stability/instability 

In order to discuss the stability of shell source in realistic stars, 
it is necessary to understand the behaviour of the thermody- 
namic quantities such as ap, ap and c* in various equation-of- 
state regimes. 

For non-relativistic partially degenerate conditions, elec- 
tron density and pressure are given by 

n e = -?— = ^(2m c kT) 3 / 2 F 1/2 ( V ) (15) 



and 



P 



pkT 2F 3/2 (7?) 
fjL a m-R 3F 1 / 2 (r]) 



(16) 



where F 3 / 2 and F x / 2 are the Fermi-Dirac functions and -q := 
ip/kT is the degeneracy parameter (e.g. Clayton 1968). The 
equation of state is thus given as: 



P = P e + P r 



pkT , pkT 2F 3/2 ( V ) f l_ aT4 



pim a /J e m H 3F 1/2 (r?) 



With the definition of := P s /P, the above equation can be 

rewritten as 

p c p, c m H 0P 



F( V ) 



Ml 



pkT 



(18) 



where F := 2F 3 / 2 /3F 1 / 2 . Differentiating this equation, we get 



1 ——dhiF=\dh\P-- — —d\nT-d\np (19) 

F + p c /pi P 

where (<91n/?/cUnP) T = (1 - 0)10 and (<91n/3/<91nT) P = 
4(/3 — l)//3 were used. From Eq. dl9l and by differentiating 
Eq. (1151 . we can obtain ap and ap as a function of 77: 



ap = 



drj 



d\nF 1/2 
dr) 



F 



dlnF 



(20) 



F + Hc/m d-q 



and 



4-3/3dlnF 1/2 



ap 



ft 



di] 



d\nF 1 / 2 



3 


F 


dlnF 


2F + 


Mc/mi 


dr] 


F 


dlnF 



(21) 



dq 



F + n e /fj,i d-q 



When j3 = 1, we get ap=ar=l for an ideal gas (?/ — » — oo, 
Fig.Q, and ap=0.6 and ap=0 for complete degeneracy (77 — > 
00). In the case of (3 7^ 0, we have ap = 1//3 and ap = 
(4 — 3/3) //3 in a non-degenerate gas (cf. Kippenhahn & Weigert 
1990), since F(rj) — » 1 and dhiF/d-q — > when r\ — > —00. 
The specific heat at constant pressure is given by cp = 
du\ p f^E.) (Kippenhahn & Weigert IT9901 . From 



\8TJP p L \u± j p 

u = (§P g + 3P r ) jp = (3 - |/3)P/a we have 



cp 



and 



V ad = 



P 

PP 

Pap 

PP^ = (| + |(1 -P))a T + 6(1-/3)- 



| + |(i-/3)]"/ -"('- ') 



ap 



(22) 



(23) 



Since P/pT in Eq. 1221 can also be given as a function of 77 
from Eq. J18I . the gravothermal specific heat c* can be obtained 
as a function of r\ for a specified 0. 

The value of can be obtained as a function of 77 at a given 
temperature from the equation of state. In Fig. ^, we show 

as a function of 77 for 4 different temperatures which may 
be relevant to the helium shell burning. /i c = 2 and p c /pi = 
0.38 have been assumed. In this temperature range, reaches 

1 when 77 > since the gas pressure becomes dominant. 

As discussed in Sections O.lt fil the values of ap and ap 
and the sign of c* affect the stability of a shell source critically. 
Fig-E> andQJ show ap and ap as a function of 77, obtained 
using the values of given in Fig.[^. The case when the radi- 
ation pressure is neglected (i.e., = 1) is also shown in these 
figures. Since changes in the total pressure become more sen- 
sitive to changes in temperature than in density as radiation 
pressure becomes significant, we find larger values of ap and 
ap when radiation pressure is considered than in the case of 
= 1. Since V a( j varies little compared to ap and ap, the 
sign of c* is mainly determined by the sign of a s ap — 1 (see 
Eq. ©). 
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Fig. 1. Thermodynamic quantities as a function of the degeneracy parameter r/ for 4 different temperatures as indicated. Here, 
T§ denotes temperature in unit of 10 8 K. As a reference, the case of j3 = 1 is also shown for ap, ap and V a d- (a) - Ratio 
of gas pressure to total pressure, 0, as a function of the degeneracy parameter, (b), (c) - ap(:— (d\n. p / dh\ P)t) and ax(:= 
(d In p/d In T)p) calculated from Eq. J20l > and Mil using j3 given in (a), (d) - V a d(:= (91nT/91nP) a d) calculated from 
Eq. d23l using j3 and «t given in (a) and (c) respectively. 



Figure|2]shows the lines of zero gravothermal specific heat 
C*, in the plane spanned by the degeneracy parameter rj and the 
relative thickness of the shell source D/r s , for 4 different tem- 
peratures. These lines have been determined using the values of 
(3, oct, Oip and V a d given in Fig.^ For (3 — 1, the thickness of 
the shell source should be larger than 0.37r s for the gravother- 
mal specific heat to be negative in the non-degenerate region 
(rj < —2). As the electron degeneracy becomes more signif- 
icant, the minimum thickness which gives c* < becomes 
larger. 

In the non-degenerate region, c* remains negative for much 
smaller shell thicknesses when the effect of the radiation pres- 
sure is considered. This is due to the fact that ap and ap 
become large when radiation pressure is significant. Since 
shell sources are generally stable when c* < as discussed 
in Sect. 13.11 this indicates that radiation pressure serves as 
a stabilizing factor for a shell source. The stabilizing ef- 
fect of radiation pressure, which has already been discussed 
by many authors (e.g. Dennis [19711 Stothers & Chin 1972 



Sackmann ll977l l. is due to the fact that a small increase in 
temperature induces considerable expansion of the shell source 
when radiation pressure is significant (cf. Eq. 10)). 

3.1. The case ofc*<0 

When c* < 0, the following two conditions are fulfilled from 
Eq. ©: 

a s ap — 1 > and a s ap — 1 < V a da s aT • 

The latter condition can be rewritten as ap — V a dCtT < l/a s - 
In a non-degenerate gas (i.e., dlnF/dr] — » 0), we have ap — 
V a dC*T — ► 0.75 from Eqs. J2CH . J2H and (J23j, while we have 
ap — VadOr — » 0.6 in a completely degenerate gas where 
ap = 0.6 and ap = 0. It follows that the inequality a s ap — 
1 < V a da s aT is always satisfied for all realistic values of a s . 
Therefore, the condition c* < reduces to a s ap — 1 > 0. With 
this condition, the term ar/(a s ap — 1) in Eq. dl4t is positive, 
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given by Harris et al. ( 1983) has been used in the calculation. 



Fig. 2. Lines of vanishing gravothermal specific heat c* in the 
plane spanned by the degeneracy parameter 77 and the relative 
thickness of the shell source D/r s , for four different temper- 
atures. These lines are computed by using the thermodynamic 
quantities given in Fig. ^ The line labeled j3 = 1 shows the 
case where radiation pressure is neglected. 



and it is generally cte > since v dominates over kt and k p 
in most cases. 

For example, let us consider a helium shell source pow- 
ered by the 3a-reaction. Figure shows v as a function of 
temperature, calculated from the 3a reaction rate given by 
Harris et al. (1983 1. The density dependence A is equal to 2 
regardless of temperature. In AGB stars or in accreting white 
dwarfs, the temperature of the helium shell source varies ac- 
cording to the core mass or the accretion rate within the range 
1 < T/10 8 K < 3, giving v = 40 •••12. Figure H shows 
the contour levels of kt and k p in the T — p plane, calcu- 
lated from the OPAL opacity table by Iglesias & Rogers ( 1996). 
The helium abundance has been assumed to be 0.66 in this 
calculation since the energy generation reaches its maximum 
around this value in a helium shell source. This figure shows 
that kt > — 1 and k p <C 1 in the given temperature and 
density range. Therefore, cte > is guaranteed as long as 
a.T/(a s ap — 1) > 0. The situation is very similar even in hy- 
drogen shell sources if they are powered by the CNO cycle. In 
conclusion, the shell source is generally stable (i.e., a < 0) un- 
der the condition c* < 0, i.e. for negative gravothermal specific 
heat. 



3.2. The case of c* >0 

For c* > 0, we have 

a s ctp — 1 < . 




Fig. 4. The logarithmic derivatives of the opacity with respect 
to to temperature and density: kt = (<91nft/<91nT) p and 
k p = (d\nK/d\np)T- The OPAL opacity table by Iglesias 
& Rogers ( 1996 1 has been used in the calculation. Y = 0.662, 
Xq = 0.286 and Z = 0.02 have been assumed. 



In this case, we have aT/(a s ap — 1) < 0. In general, a shell 
source is prone to thermal instability when aT/(a s ap — 1) < 
since any density decrease due to expansion will lead to a fur- 
ther temperature increase (see Eq. @). However, if the abso- 
lute value of aT/(a s ap — 1) is large enough, a small increase 
in temperature induces such a large decrease in density that the 
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nuclear energy generation may be considerably reduced. In this 
case, we have cte < and the shell nuclear burning is thermally 
stable. 

This can be realized when a^ becomes very large or/and 
when \a s ap — 1| <C 1. A large ax can be obtained if radia- 
tion pressure becomes substantial (Fig.^), showing again its 
stabilizing effect. When |a s ap — 1| -C 1 is satisfied, we have 
Sp/p — —SP/P ~ apSP/P. This means that the shell source 
behaves as if it were isothermal when a s ap ~ 1. As a conse- 
quence, the influence of the temperature perturbation becomes 
negligible compared to that of the density perturbation, result- 
ing in cte < 0. For an ideal gas with f3 = 1, this condition 
means a s — V. the shell source is stable if its thickness D is 
not much less than 0.37r s , even when c* > (see Fig.|5}- This 
mode of stability has been already noted by Fujimoto ( 1982b i: 
he found that a very large c* can render a shell source stable, 
which is the case when \a s ap — 1| -C 1 (Eq. (|6|l). As shown in 
Sect. |5] this situation may occur in the helium shell sources of 
AGB stars and accreting white dwarfs. 



3.3. Infinitely thin shell sources 



o, 



cp 



(24) 



If a shell source is infinitely thin, we have a s 
and 

a = v — 4 + «t — «t(A + K p ) 

The shell burning can be stable if the dependence of the nuclear 
reaction on temperature is weak or/and when radiation pressure 
is significant, in which case ar becomes large (cf. Fig.|5}. An 
astrophysical application for such a case is shell burning in ac- 
creting neutron stars, for which the compactness of the grav- 
itating body enforces a constant pressure in the shell source, 
such that a s — > 0. In fact, the criterion for thermal instability in 



this case is very similar to that given by Fujimoto et al. ( 198 1 1, 
who studied shell flashes on accreting neutron stars. Numerical 
calculations by Taam et al. ( 1996 1 show that the helium shell 
burning in a neutron star can be stable when the temperature 
of the helium shell source is as high as 5.6 x 10 8 K, for which 
the temperature sensitivity of the 3a reaction is considerably 
weakened (v ~ 5, Fig. [3}. This can be easily understood by 
our stability criterion with Eq. J24l > (see also Fig.Qand discus- 
sions given in Sect.|4j. 
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Fig. 5. Loci of ct and c* for their sign in rj — D/r s plane when 
T = 1.7xl0 8 K. The solid line separates the thermally unstable 
(ct > 0) from the stable (ct < 0) region, while the dashed line 
separates the region with negative c* from that with positive c*. 
The same chemical composition as in Fig.|4]is assumed. 
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3.4. Complete degeneracy 

In completely degenerate gas (i] 
ct becomes 



7] 



oo), we have ar = and 



Kt- 



(25) 



Since heat in a strongly degenerate gas is transfered mainly 
by electron conduction, with kt — 2 (e.g. Kippenhahn & 
Weigert 1990), any nuclear reaction with v > 2 will give 
ct > 0, resulting in unstable nuclear burning. An example 
for the shell burning in strongly degenerate conditions is pro- 
vided by the sub-Chandrasekhar mass progenitor model for 
Type la supernovae. In this model, helium ignites at the bot- 
tom of a helium layer of some tenth of solar mass accumulated 
on a CO white dwarf with low accretion rates (e.g. Woosley 
& Weaver 1994 1. At the ignition point, the helium envelope is 



Fig. 6. Loci for the sign of ct at 4 different temperatures: T 
1.2, 1.7, 2.2 and 3.0 x 10 8 K in <q - D/r s pla— 



lane. 



highly degenerate (77 > 10, e.g. Yoon & Langer 2004 1, render- 
ing complete degeneracy as a good approximation. 

4. Instability diagrams 

Figure |5] gives an example of the stability criterion for a He 
shell source in 77 — D/r s parameter space, using a specific 
chemical composition as indicated in the figure caption, which 
are typical for helium shell sources in AGB stars or in accreting 
white dwarfs. 
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Fig. 7. Stability conditions for a helium shell source in the 
density-temperature plane. The solid line separates the ther- 
mally unstable (a > 0, lower right part) from the stable (a < 0) 
region, for 5 different relative shell source thicknesses (i.e., 
D/r s = 0.0, 0.1, 0.2, 0.3 and 0.4). The dotted contour lines de- 
note the degeneracy parameter r\. X^ B = 0.662, Xq — 0.280 
and Z = 0.02 have been assumed. 



This figure shows the two modes of stability as discussed in 
sections IXTl and l3~2l one with c* < & cte > and the other 
with c* > & cte < 0. As expected, a shell source becomes 
more susceptible to thermal instability when it is thinner and 
more degenerate. For r\ > 2, where the electron degeneracy 
is significant (F(rj) > 1.5), the shell source can not be stable 
even with D/r s > 0.5, which is because a pressure change 
becomes less sensitive to a temperature change as the electrons 
become more degenerate. On the other hand, the shell source 
can be stable even with D/r s = when r\ < —5.8 (see also 
Fig. [7J, where the radiation pressure is dominant over the gas 
pressure. 

As shown in Fig.|4] kt and k p remain to be small compared 
to v in the given temperature and density range. The quantities 
Qp and ar do not depend strongly on /j c /mi- The main factor 
which determines the stability at given 77 and D/r s is thus the 
temperature of the shell source. Figure |6] shows that the sta- 
bility criterion is more relaxed for higher temperature. This is 
because the temperature dependence of the energy generation 
is weakened (Fig.|5J, and because radiation pressure becomes 
more significant at a given degree of degeneracy as the temper- 
ature increases. This temperature dependence of the stability 
of shell burning explains, why the accretion rates which allow 
steady burning solutions depend on the mass of accreting white 
dwarf (e.g. Fuiimoto fT982bl Nomoto & Kondo ll99ll Cassisi et 
al.|1998). I.e., higher accretion rates are required to have stable 
shell burning for more massive white dwarfs, in which shell 
sources are thinner. This is mainly because shell sources be- 
comes hotter with higher accretion rates if the accreted matter 
is burned in a steady way. 




2 3 

log p [g/cm 3 ] 

Fig. 8. Stability conditions for hydrogen shell source pow- 
ered by the CNO cycle, in density and temperature plane. The 
solid line separates the thermally unstable {a > 0) from sta- 
ble (a < 0) region at 5 different relative thicknesses of shell 
source (i.e., D/r s = 0.0, 0.1, 0.2, 0.25 and 0.3). A shell source 
with each D/r s is expected to be stable in the upper region of 
the each locus. The dotted contour lines denote the degeneracy 
parameter 77. X K = 0.502, A" Hc = 0.478 and Z = 0.02 have 
been assumed. Note that the hot CNO cycle is not considered 
in the calculation. 



In Fig. we show the stability conditions in the density- 
temperature plane for 5 different relative shell source thick- 
nesses. As discussed in Sect. 13.31 this figure shows that a shell 
source can be stable even when it is infinitely thin (D/r s = 
0), if the temperature is high enough, since with high temper- 
ature the dependence of the energy generation on temperature 
becomes weakened and the role of radiation pressure becomes 
important, giving a large olt- The stability conditions are re- 
laxed as D/r s increases. The conditions for stability in hy- 
drogen shell sources show a similar behaviour (Fig. [8}, where 
the shell source is assumed to be powered by the CNO cycle. 
I.e., thinner, cooler and more degenerate shell sources are more 
prone to thermal instability. Note, however, that in the calcula- 
tion, we did not consider the so called hot CNO cycle, which 
may be active at high temperatures (T > 8 x 10 7 K). 

5. Applications 

In this section, we apply the stability criterion to various nu- 
merical stellar models: a 3.0 M Q AGB star, helium accret- 
ing CO white dwarfs, a cooling helium white dwarf with a 
thin hydrogen envelope, and a 1 .0 M Q giant. All models have 
been computed with a stellar evolution code (Langer 1998 and 
references therein) which solves the hydrodynamic form of 
the stellar structure equations (Kippenhahn & Weigert 1990 1. 
Opacities are taken from Iglesias & Rogers ( 1996 1. We define 
the thickness of a shell source D such that the energy genera- 
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Table 1. Properties of the computed helium accreting white 
dwarf models. The first column denotes the model sequence 
number. The second column shows the initial mass and the 
third gives the accretion rate. The forth column gives the 
amount of accreted material until the onset of the helium shell 
instability. The value of a estimated in the shell source of the 
last analyzed model is given in the last column. 
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0.0212 


-0.429 


No.2 


0.800 


5.0- 10~ 7 


0.0238 


-0.291 


No. 3 


0.800 


7.0 ■ 10~ 7 


0.1141 


-0.598 


No.4 


0.800 


8.0 ■ 10~ 7 


0.1455 


-0.471 


No.5 


0.998 


8.0 ■ 10" 7 


0.0028 


-0.114 


No.6 


0.998 


1.0 -lO -6 


0.0048 


-0.035 



4.300 4.305 
time (10 8 yrs) 



1.80 | ' ' 
1.50 r 
1.40 r 
1.30 j- 

1.20 r 
4.295 



(d) 



(c) 



0.38 
0.36 
0.34 
0.32 

0.30 



4.300 4.305 4.310 
time (10 8 yrs) 



-1.2 -1.0 -0.8 -0.6 

n 



Fig. 9. Evolution of the helium shell source in a 3 M Q AGB star, 
(a) Nuclear luminosity due to helium burning as a function of 
time, for shortly before and during the onset of thermal pulses. 
The mean values of c*, a, and temperature in the shell sources 
of the models marked by a filled circle are shown in (b), (c) 
and (d) respectively, (e) Evolution of the shell source in the r\ — 
D/r s plane. The dotted line corresponds to a = computed for 
the physical conditions in the shell source of the last analyzed 
model. 



tion rate at each boundary has a certain fraction (:= /b) times 
its peak value. Values of D/r s in stellar models would be af- 
fected by different choices of /b- For instance, /b = 1CP 3 re- 
sults in a larger D/r s by about 15 % compared to the case of 
/b = 1CP 2 , in helium accreting white dwarf models presented 
in Sect. 15.21 In our applications, we choose /b = 2 x 10~ 3 be- 
cause this particular choice gives the best prediction for the on- 
set of thermal instability of shell burning for all different types 
of stellar models discussed below. The mean values weighted 
by the energy generation rate over the shell source are used 
for density, temperature and chemical abundances to estimate 
physical quantities which appear in Eq. d!2i . 

5. 1. Helium shell burning on the AGB 

It is well known that asymptotic giant branch (AGB) stars suf- 
fer thermally unstable helium shell burning in their advanced 
evolution (e.g. Iben & Renzini 1983). We investigate this phase 
by computing a 3 Mq evolutionary model with Z = 0.02, start- 
ing at the zero age main sequence (Langer et al. 119991 . After 



passing through hydrogen and helium core burning, the star en- 
ters the helium shell burning stage. Fig. |9^ shows the nuclear 
luminosity due to helium burning as a function of time at the 
onset of the thermal pulses (T ~ 4.306 x 10 8 yrs). We an- 
alyzed those models indicated by a filled circle in Fig. [9^, in 
which the helium shell burning is still stable, and computed c* 
and a for each model. (Fig.[9j> and Fig.|9j;). During the stable 
shell burning stage, the mean temperature in the shell source is 
nearly constant at T ~ 1.4 x 10 8 K (Fig.|9jl). The shell source 
becomes thinner and more degenerate as the CO core mass in- 
creases as shown in Fig. [9^. 

The gravothermal specific heat of the shell source is posi- 
tive in the chosen models, indicating that the shell burning is 
kept to be stable, until the onset of the thermal pulses by the 
second mode of the stability discussed in Sect. 13.21 However, 
it becomes smaller as the CO core mass increases, mainly be- 
cause a s ap-l deviates from 1 more significantly as D/r s de- 
creases (see Eq. (|6j), and partly because degeneracy in the shell 
source becomes stronger (Fig. |9^). The thermal pulses start 
when the relative thickness of the shell source becomes as low 
as 0.316. 

5.2. Helium shell burning in accreting white dwarfs 

Here, we consider the helium shell source in hot helium accret- 
ing CO white dwarfs. Two different initial white dwarf models 
are considered: M; n ;t = 0.8 M Q with logi s /L Q = 2.948 and 
Mi n it 0.998 M Q and logL s /L Q = 3.379, where M init and L s 
are the initial mass and surface luminosity, respectively. The 
accreted matter is assumed to be helium enriched such that 
Y = 0.98. We have performed simulations with 6 different 
combinations of initial mass and accretion rate (see Table 1). 
The helium shell burning is initially stable in each model se- 
quence, and eventually turns unstable as the white dwarf mass 
increases (see Fig.[K)]for an example). The value of o for the 
shell source of the last analyzed model in each sequence is 
given in the last column of the table. In all sequences, the last 
analyzed models, from which the thermal pulses grow signif- 
icantly, have \a\ < 0.6, indicating that our stability criterion 
predicts the onset of the instability with good accuracy. 
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Fig. 10. Evolution of the helium shell source in an accreting 
white dwarf with an initial mass of 0.998 M© and a constant 
accretion rate of 1CP 6 Mgyr -1 . (a) Nuclear luminosity due 
to helium burning as a function of the total white dwarf mass. 
The mean values of c* , a, and temperature in the shell sources 
of the models marked by a filled circle are shown in (b), (c) 
and (d) respectively, (e) Evolution of the shell source in the r\ — 
D/r s plane. The dotted line corresponds to a = computed for 
the physical conditions in the shell source of the last analyzed 
model. 



As an example, results for sequence No. 6 are shown in 
Fig-El In Fig.llOh. the nuclear luminosity due to helium burn- 
ing is given as a function of the total white dwarf mass, which 
serves as a linear measure of time. Thermal pulses are induced 
when the white dwarf mass reaches ~ 1.028 M Q . The evolu- 
tion of c* and a shown in Fig.llOb andllOt is very similar to that 
in the AGB model discussed above. The mean temperature in 
the shell source is about 2 x 10 8 K (Fig. HOB ). The shell source 
becomes thinner and more dense as the white dwarf mass in- 
creases (Fig. HOt . The last analyzed model is at the edge of 
stability as shown in Fig.llOt and the thermal pulses grow sig- 
nificantly from this point. 
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Fig. 11. (a) Nuclear luminosity due to hydrogen burning at the 
onset of the hydrogen shell flash in the hydrogen envelope in 
a cooling helium white dwarf of 0.324M Q . (b) Contour levels 
of the perturbation growth time scale in the rj — D/r s plane, 
obtained from Eq. J13I . Here, a and r t h have been estimated 
for the physical conditions in the shell source of the marked 
model. The solid line corresponds to a = 0. 



5.3. Hydrogen shell burning in cooling helium white 
dwarfs 

Here, we discuss a hydrogen shell flash in a cooling helium 
white dwarf model (cf., Driebe et al. I1999> . which has been 
obtained by evolving a 2.4 M Q zero age main sequence star 
with Z = 0.02 in a binary system with the initial orbital pe- 
riod of 1.76 days. Its companion was a 2 M Q zero age main 
sequence star. It loses most of the hydrogen-rich envelope dur- 
ing the hydrogen shell burning stage through so-called case B 
mass transfer, ending in a helium white dwarf of M — 0.324 
M with a hydrogen envelope of M ~ 0.005 M© when the 
mass transfer stops. The luminosity and surface temperature at 
this moment are 158.5 L© and 1.4 x 10 4 K, respectively. 

The hydrogen shell burning is weakened progressively as 
the white dwarf cools (Fig. [Tiki, At t ~ 5.905 x 10 8 yr, 
further weakening of the shell burning stops due to the in- 
crease in density. The hydrogen burning increases slowly as 
the helium core contracts further and ends in a strong flash 
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at t ~ 5.98 x 10 8 yr. The thickness and degree of degen- 
eracy of the shell source have been computed for the model 
marked by a filled circle in Fig.lllh. and positioned in D/r s — r\ 
plane in Fig. II lb . The corresponding temperature and density 
of the shell source are 2.13 x 10 7 K and 6.91 x 10 2 g cm" 3 , 
respectively. The shell source at this point is in the unstable 
regime as expected from the fact that the nuclear luminosity 
reaches its peak very soon thereafter. The contour levels for 
the perturbation growth timescale from Eq. Jl 3I > are also given 
in the same figure. The perturbation growth time estimated at 
the point is about 2.4 x 10 4 yrs, which is in good agreement 
with the time interval from the marked model to the peak in 
I H (^3x 10 4 yrs). 

5.4. Hydrogen shell burning in a 1.0 M giant 

It is a matter of debate whether shell flashes in hydrogen shell 
sources of low-mass sub-giants and giants could be a possible 
cause for anomalies of chemical abundances in stars of glob- 
ular clusters (e.g. Denissenkov & van den Berg 2003 1. Stellar 
evolution models indicate that hydrogen shell sources in sub- 
giants or red giants are usually thermally stable. Here, we ana- 
lyzed the physical conditions of the hydrogen shell source in a 
1.0 Mq model with solar metalicity. 

Fig.ll2h shows the evolution of the hydrogen shell source in 
the density-temperature plane, from the beginning of the sub- 
giant phase until the core helium exhaustion. Note that, from 
the beginning of the sub-giant until the core helium flash, the 
shell source becomes hotter and less dense (Fig.ll2h). while the 
relative thickness of the shell source remains larger than 0.3 
(Fig. 112b ). Once the core helium flash occurs, both the tem- 
perature and density of the shell source decrease rapidly due 
to the expansion. Throughout the considered stages, the evolu- 
tionary track of the shell source remains in the thermally stable 
regime (cf. Fig. 112b). This result indicates that the onset of the 
thermal instability in hydrogen shell sources of low-mass sub- 
giants and giants is unlikely. 
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Fig. 12. Evolution of a 1.0 star with solar metalicity from the 
beginning of the sub-giant phase until the exhaustion of core 
helium burning, (a) Evolution of the hydrogen shell source in 
the density - temperature plane. The thick solid lines corre- 
spond to a = for D/r s = 0.3 and 0.4 (cf. Fig.|SJl. The evo- 
lutionary stages at different epochs are indicated with short ar- 
rows. The dotted contour lines give the degeneracy parameter 
r\. (b) The relative thickness of the shell source as a function 
of time, (c) The value of a in the shell source as a function of 
time. 



6. Conclusion 

We have developed a stability criterion for thermonuclear shell 
sources which can be easily applied to numerical stellar models 
in a quantitative way. The stabilizing/destabilizing factors have 
been identified unambiguously: a shell source is less prone to 
the thermal instability when it is thicker, less degenerate and 
hotter (Sect.|5J Figs.|6j0and|8|l. If a shell source is sufficiently 
thick at a given degree of degeneracy, the gravothermal spe- 
cific heat becomes negative and any excessive energy gain is 
consumed mostly for the expansion work, preventing thermal 
runaway. It is shown that the shell source can remain stable 
even with a positive gravothermal specific heat as long as tem- 
perature at a given degree of degeneracy is high enough such 
that the sensitivity of nuclear reactions to temperature is signif- 
icantly weakened or/and that radiation pressure becomes sub- 
stantial. In such a case, the energy loss rate by radiation be- 
gins to dominate over the additional nuclear energy production 
rate. Interestingly, this effect allows even an infinitely thin shell 
source to be thermally stable (Sect. l3~3l Fig.0and[8}, which 



may have interesting consequences in accreting neutron stars 
(Sect 1331 cf. in't Zand et al. 120051 . 

We have applied the stability criterion developed in this 
study to a 3 M AGB star, to helium accreting CO white 
dwarfs, to a helium white dwarf with a thin hydrogen shell, 
and to a 1.0 M giant (Sect.|5Jl. In spite of the simple approach 
of using homology assumption in the shell source and not rely- 
ing on the full linear stability analysis, it is demonstrated that 
this criterion can predict the onset of thermal pulses or flashes 
with a reasonably good accuracy. 

In conclusion, being simple and robust, our stability cri- 
terion may serve as a useful tool for analyzing shell burning 
in various types of stars. One example can be found in Yoon 
et al. (2004 1, who discuss the effects of rotationally induced 
chemical mixing on the behaviour of helium shell burning in 
accreting white dwarfs. 
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